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We propose a technique, which produces nearly complete ionization of the population of a dis- 
crete state, coupled to a continuum by a two-photon transition via a lossy intermediate state, whose 
lifetime is much shorter than the interaction duration. We show that using counterintuitively or- 
dered pulses, as in stimulated Raman adiabatic passage (STIRAP), wherein the pulse coupling the 
intermediate state to the continuum precedes and partly overlaps the pulse coupling the initial and 
intermediate states, greatly increases the ionization signal and strongly reduces the population loss 
due to spontaneous emission through the lossy state. For strong spontaneous emission from that 
state, however, the ionization is never complete because the dark state required for STIRAP does 
not exist. We demonstrate that this drawback can be eliminated almost completely by creating a 
laser-induced continuum structure (LICS) by embedding a third discrete state into the continuum 
with a third, control laser. This LICS introduces some coherence into the continuum, which enables 
a STIRAP-like population transfer into the continuum. A highly accurate analytic description is 
developed and numerical results are presented for Gaussian pulse shapes. 
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PACS numbers: 32.80.Qk, 32.80.Fb, 32.80.Rm, 33.80.Rv 



I. INTRODUCTION 

Stimulated Raman adiabatic passage (STIRAP) is a 
simple, robust and efficient technique for complete popu- 
lation transfer (CPT) in three-state quantum systems 
In this technique, the population is transferred adiabati- 
cally in a Raman transition, from an initially populated 
state ipi via an intermediate state ip2 to a target state ip3 
by two pulsed fields, pump and Stokes, whose frequencies 
are maintained on two-photon resonance between states 
ipi and ip3- If the pulses are ordered counterintuitively, 
the Stokes before the pump, then the dark state is associ- 
ated with state ipi initially and state ^3 in the end, thus 
providing an adiabatic route from ip\ to ip3. A unique and 
remarkable feature of STIRAP is that during the trans- 
fer the population remains trapped in a dark state, which 
is a time-dependent coherent superposition of states ipi 
and ip3 only and does not involve the intermediate state 
ip2- State ip2 therefore remains unpopulated during the 
transfer and its properties, including possible population 
decay, are largely irrelevant for STIRAP. A very large 
detuning or loss rate 0, however, do affect STIRAP 
and reduce its transfer efficiency. 

The simplicity, efficiency, and robustness of STIRAP 
have attracted much attention, which has resulted in nu- 
merous applications in a variety of quantum systems. 
Among them, we mention population transfer via a con- 
tinuum, wherein the discrete intermediate state tp2 is re- 
placed by an ionization continuum [fa-|6[ • Because a dark 
state that links adiabatically states i[>i and tp3 still forms, 
high transfer efficiency is still possible. However, various 
specific features of the continuum that make it substan- 
tially different from a discrete state (or a manifold of 
discrete states) reduce the transfer efficiency. Several sce- 



narios have been proposed to reduce the effects of these 
continuum features, resulting in increase of the transfer 
efficiency well above 50%. STIRAP via continuum has 
recently been demonstrated experimentally Q. 

The extension of STIRAP to include an initial or final 
continuum was at first investigated for the photoasso- 
ciation process producing cold molecules starting from 
laser-cooled atoms 0. Photoassociation is based on the 
laser-driven transitions of a continuum-bound-bound sys- 
tem. Vardi et al 9] have extended the STIRAP technique 
to the cases of initial or final continuum, proving that 
transfer with a good efficiency is possible. The quantum 
transfer in a three-level system based on bound-bound- 
continuum transitions arises also in the ionization of a 
Bose- Einstein condensate (BEC). Ionization of a rubid- 
ium BEC from the ground state through a two-photon 
ionization (TPI) scheme was explored experimentally [lOj 
and theoretically fill ]. 

In the present paper, we follow this latter idea and 
introduce a continuum in STIRAP by replacing the fi- 
nal state ip3 by a continuum. The objective, compared 
to STIRAP via a continuum, is also changed drastically: 
instead of trying to avoid ionization, here we aim at max- 
imizing ionization. As in STIRAP we also aim at mini- 
mizing the transient population of the intermediate state 
ip2, in order to avoid (possibly strong) decay to other 
discrete states via spontaneous emission or unwanted ex- 
citation to other states. States ipi and 4>2 are linked by 
a pump pulse, and state tp2 is connected to the contin- 
uum by an ionizing pulse, with the ionizing pulse arriv- 
ing before the pump. The challenge here is that, unlike 
STIRAP via continuum, a dark state cannot be formed 
between the initial state ipi and the continuum states, 
because a flat continuum is an incoherent medium. 



To this end, we propose to use a laser-induced contin- 
uum structure (LICS) [7| created in the continuum by 
embedding an ancillary, control state ip c into the contin- 
uum by a third, control laser. This LICS creates some 
coherence into the continuum, which makes it possible to 
create a quasi-dark state, and thence a STIRAP-like pro- 
cess into the continuum. This LICS-STIRAP allows us to 
produce almost complete ionization, with negligibly small 
population losses from state 02, even when the ip2 life- 
time is much shorter than the interaction duration. Ex- 
perimental verification of this LICS-STIRAP technique 
will open opportunities for many applications, for in- 
stance efficient photoionization of a BEC without atomic 
excitation in intermediate atomic states that will improve 
the ion production obtained in the experiment of ref. 
fiol j . In addition, an increase in the ionization efficiency 
for cold atoms will improve the brightness in magneto- 
optical-trap based sources producing either electron [l2[ 
or ion [13J beams. 

The structure of this paper is as follows. In Sec. |TT] 
we define the problem and show that a counterintuitive 
pulse order of the pump and ionizing pulses suppresses 
fluorescence and is favourable for ionization, even with- 
out a control laser. In Sec. IIIII we add LICS to the 
scenario and derive a very accurate analytic approxima- 
tion to describe this LICS-STIRAP technique. Section 
IIVI provides illustrations of the proposed ionization tech- 
nique and comparison of the analytical and numerical 
results. Section [V] presents a summary. 



II. STIRAP INTO CONTINUUM 

We first consider TPI of an atom, initially in state -01, 
coupled to the ionization continuum via state 02 > as il- 
lustrated in Fig. [1] The transition 1-2 is driven by a 
pump laser pulse with Rabi frequency fl(t) and detuning 
A, and state ip2 is connected to the continuum by a sec- 
ond laser pulse with a time-dependent rate Fi(i). State 
02 can decay irreversibly via spontaneous emission (or 
other mechanisms) to other states with a constant rate 
T; we shall refer to the respective signal as the fluores- 
cence signal, 



F = 



ionization 
laser 




pump 
laser 



FIG. 1: (Color online) TPI scheme. An initially populated 
discrete state 0i is coupled to another discrete state 1^2 by 
a resonant, or nearly resonant, pump laser field. State 02 
is coupled to the ionization continuum by a second, ionizing 
laser pulse. State 02 can decay irreversibly to other states. 



02. This objective reminds one of STIRAP. The signif- 
icant difference here is that state 03 is replaced by a 
continuum. 

The equation that describes the dynamics of the sys- 
tem (in units h = 1) is the Schrodinger equation, 



ij t <t) = H(i)c(t). 



(3) 



Here c(t) — [ci(t), c<i{t)\ is the column- vector with 
the probability amplitudes C\(t) and C2(t) of states tpi 
and 02 1 and H(i) is the Hamiltonian, obtained by adia- 
batic elimination of the continuum states and within the 
rotating- wave approximation (RWA) [7] , 



TP 2 (t)dt. 



(1) 



H(t) = \ 



25i(t) 



2A- 



fi(t) 
■2S 2 {t)-iT l {t) 



iT 



(4) 



The ionization signal is 



7 = l-P 1 (oo)-P 2 (oo)-F, 



(2) 



where P n (t) are the populations of the discrete states 
(n = 1,2). Our objective is to set up the laser pulses 
such that the ionization signal I is maximized, while the 
fluorescence signal F is minimal. 

We wish to design a recipe to maximize ionization 
when the loss rate T is large compared to the interaction 
duration T. In other words, we wish to ionize the atom, 
without exciting it, despite being on resonance with state 



where Si(t) and 62(4) are the Stark shifts of states ipi 
and %p2, produced by virtual excitation to other atomic 
states. 

We are interested in situations when the ionization sig- 
nal I is large, i.e. in the non-perturbative regime. This 
implies large peak Rabi frequency and large ioniza- 
tion rate Ti(t). Because we also assume that the loss rate 
T is fixed and large, T ^> 1/T, state ip 2 is subjected to 
strong population decay, due to both spontaneous emis- 
sion and ionization. This implies that it receives very lit- 
tle transient population and can therefore be eliminated 
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FIG. 2: Contour plots of the ionization signal 7 as a function 
of the delay of the ionizing pulse Tj and the pump detuning 
A = 0. The pulses have Gaussian shapes, Eqs. (|22a|) . The 
decay rate from state tp2 is T = 100/T, the peak ionization 
rate is Tio — 50/T, the peak Rabi frequency is fin = 50/T, 
and the Stark shift is S = 0. 



FIG. 3: Contour plots of the ionization signal I as a function 
of the delay of the ionizing pulse Ti and the peak ionization 
rate Fm. The pulses have Gaussian shapes, Eqs. (|22a|) . The 
pump detuning is A = 0, the decay rate from state ip2 is 
r = 100/T, the peak Rabi frequency is fin = 50/T, and the 
Stark shift is S = 0. 



adiabatically. Hence we find after simple algebra 



Pi(t) « exp 



[r l (t') + r] 2 + 4[A + ^(t')] 



W m) 

[Ti(t) + T} 2 + 4[A + S(t)f W ' 



(5b) 



with S(t) = Sa(t) — Si(t). Now we use these formulas to 
examine the possibilities of how to minimize fluorescence 
F, i.e. Pi(t), and simultaneously maximize ionization /. 

The first choice is to use a large detuning A. In- 
deed, this will reduce the population of state ip2, because 
Pi(t) ~ A~ 2 for large A [see Eq. (|5bj) ]: however, this 
decrease will be accompanied by an increase in the pop- 
ulation of state ipi, see Eq. (|5ap . In result, the increase in 
the ionization will be little, if any. A similar conclusion 
applies to the Stark shift S(t) , which can be induced by 
the ionizing laser, or, if needed, by an additional far-off- 
resonance laser, as in SCRAP technique flij ]. 

A straightforward alternative is to increase the magni- 
tude of the ionizing pulse Ti(t) alone. Then the popula- 
tion ([FBI of state ifi? decreases as T~ 2 , as vs A; however, 
the increase in the population (|5aj) of state ip\ is smaller 
vs Ti than vs A. In result, the ionization signal will 
increase more markedly when increasing Ti . 

A closer inspection of Eqs. ([3]) suggests that one can 
decrease P2(t), without increasing P\(t) (implying thence 
a net increase of the ionization I) , by delaying the pump 
pulse Cl(t) with respect to the ionizing pulse Ti(t). In- 
deed, it is obvious that the pump pulse Q(i) must not 
arrive before the ionizing pulse T(t), because then the flu- 
orescence will deplete the population even before ioniza- 
tion has the chance to begin; mathematically, this implies 
large values for the fractions in Eqs. ([5]), with resulting 



small population of state ipi and large fluorescence sig- 
nal F. In contrast, if the pump pulse f2(t) arrives simul- 
taneously, or after the ionizing pulse Ti(t), with some 
overlap, then fluorescence can only begin simultaneously 
with ionization. Moreover, if during the ionization the 
ratio H{t)/Ti(t) is very small, while the ratio n(t) 2 /Ti(t) 
is moderately large compared to 1/T, then both P\(t) 
and Pi{t) will remain small, as easily seen from Eqs. ([5]) 
when A = S = 0. Therefore, our objective of produc- 
ing ionization without excitation requires that during the 
population depletion we have 



n(t) 2 T> [Ti(t)+T] > n(t). 



(6) 



These conditions require large pulse areas over the inter- 
action duration, 

/OO /"OO 
n(t)dt^>i, / Ti{t)dt^i. (7) 
-co J — OO 

Conditions ^ suggest that the pump pulse Rabi fre- 
quency Q,(t) should be small in comparison with the ion- 
izing rate Ti(t). This can be naturally achieved, indeed, 
if the pump pulse is delayed to, but overlapped with, 
the ionizing pulse. Then ionization will occur during the 
rising edge of the pump pulse. 

These conclusions are illustrated in Figs. [5]and[3]where 
the ionization signal is plotted as a function of, respec- 
tively, the ionizing pulse delay and the pump pulse de- 
tuning A, and the ionizing pulse delay and the ionizing 
pulse intensity. These figures clearly demonstrate that 
the counterintuitive pulse order - ionizing pulse before 
pump pulse - is favorable for ionization. Figure [5] demon- 
strates also that the detuning A is of little help in respect 
to ionization, as predicted by Eqs. (JS|). 
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FIG. 4: (Color online) The scheme for ionization by LICS- 
STIRAP, which extends the scheme in Fig. fT] The additional 
discrete state %p c is embedded into the continuum by a control 
laser field, which creates a LICS in the continuum. 



As follows from the above analysis, counterintuitive 
pulse order increases ionization and suppresses excita- 
tion. However, very strong ionizing laser is needed to en- 
sure the conditions (JSJ) ■ In the following section we shall 
show that a LICS in the continuum can help ionization 
and make this process very similar to STIRAP because 
of the creation of a quasi-dark state between state ip\ and 
the LICS, with ensuing nearly complete ionization with 
moderate laser resources. 



in population transfer via continuum 4]. The quanti- 
ties F 2 = £ Q T« and T c = Y, a T c (<* = P»*,c) are 
the total ionization rates of states ip2 and ip c , respec- 
tively, which are given by sums of ionization rates in- 
duced by the pump (p), ionization (i) and control (c) 
pulses, whereas S n = J2 a S" (n — 1,2, c) are the dy- 
namic Stark shifts. The ionization widths and the Stark 
shifts are proportional to the pulse intensities I p (t), Ii(t) 
and I c (t), r°(i) = r»o/ a (t) and S^(t) = S% I a (t), where 
the parameters T" and S" depend on the particular 
atomic states and the laser frequencies. 

For simplicity, we will assume that the ionization of 
state ip2 occurs due to the action of the ionizing laser p 
only: r 2 = T 2 ; this condition can be satisfied by selecting 
appropriate atomic levels and laser frequencies such that 
the pump and control lasers do not cause direct ionization 
from state ifa. We will also assume that the ionization 
rate of state ifj c is induced by the control laser I c only: 
r c = rjj. Given the preceding assumption for r 2 it is 
clear that the ionizing laser Ii will ionize also state ip c 
(r* 7^ 0). However, because from state i/j c the ionizing 
laser points deeply into the continuum this ionization rate 
is small; moreover, as state %p c remains largely unpopu- 
lated, this ionization channel (which is actually favorable 
for our goal of maximizing ionization) does not alter the 
dynamics markedly. 

In addition, we also neglect the Stark shifts, S\ — S2 = 
S c = 0; these are important in population transfer via 
continuum [H, but do not have much effect here. 

The fluorescence signal F from state ip2 is given by Eq. 
([1]), and the ionization signal by 



I = 1 - Pi (00) - P 2 (oo) - P c (oo) - P. 



(9) 



III. STIRAP INTO LICS 

A. The system 

Let us assume now that, in addition to the scheme in 
Fig. [TJ an additional laser pulse couples a third discrete 
state ijj c with the continuum, as shown in Fig. 2J The 
dynamics of the system is again described by Eq. ([3]), 
with c(t) = [ci(t), C2(t), c c (t)] T . The Hamiltonian after 
adiabatic elimination of the continuum states and within 
the RWA reads 



2S*i fi 
2A + 2S 2 - iTj- iT 
-(q + i)VT^T c 







-(<? + i)\/TiT c 
25 + 2S C - iV c 



(8) 

where explicit time dependences are omitted for brevity. 

The constant q, called the Fano parameter [iEj . is an 
important feature of LICS. It is responsible for the asym- 
metric dependence of the ionization signal on the two- 
photon detuning Q, and it also plays an important role 



The optimal pulse order of the three pulsed fields is 
determined by the objective to maximize ionization. We 
have already come to the conclusion that the pump and 
ionizing pulses must arrive in a counterintuitive order: 
ionizing before pump. The optimal timing of the control 
pulse can be deduced from the following arguments. 

Because the objective is ionization, we must avoid pop- 
ulation transfer via the continuum into state tp c , which 
will occur if the control pulse precedes the ionizing pulse 
[H ; the control should therefore be applied after the ion- 
izing pulse. The timing of the control pulse with respect 
to the pump pulse is not so significant but these pulses 
should not be separated too much because it is obvious 
that, for LICS to have any effect, the control pulse must 
overlap significantly with the pump and ionizing pulses. 
We therefore conclude that for maximal ionization, the 
pulses should be applied in the order ionizing- control- 
pump, with a sufficient overlap between them. We shall 
therefore assume this pulse ordering in the analytical de- 
scription in the next section, which will be confirmed as 
optimal also by numerical simulations in Sec. IIVI 
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B. Analytical description 

It is appropriate to describe the evolution of the system 
in the basis of the instantaneous eigenstates of the Hamil- 
tonian ([5]). Because this Hamiltonian is non-Hermitian 
its eigenvalues e a (t) are complex valued, and the right 
eigenvectors f a (t) differ from Hermitian conjugates of 
the left eigenvectors 16]. The right eigenvectors, and the 
eigenvalues, are defined by the equation (a = +, 0, — ) 



H(t)ip a (t) = e a (t)ip a (t), 



(10) 



where ip a (t) = [fi a (i), f2a{t), h a {t)\ 
we expand the state vector as 

*(t) =b+(t)ip + (t) + b (t)ip (t) + 



. Using such states 
&-(%-(*)• (11) 



The probability amplitudes in the original basis and the 
adiabatic basis are related through the transformation 



c(i) - R(t)b(i), 



(12) 

T 



where the column vector h(t) — [b+(t),bo(t),b-(t) 
comprises the probability amplitudes of the adiabatic 
states. The columns of R are the components of ip a (t), 
( R W)„ Q = /na(*), with n = 1,2, 3 and a = +, 0, -. 
The Schrodinger equation for the vector b reads 



i-h(t) = H ad (t)h(t), 



(13) 



where H ad (t) — H a (t) 
nal Hamiltonian 



H raa (i), with an adiabatic diago- 



H°(t) = R- x (i)H(t)R(t) 



+ (*) 








e (t) 







e _(t) 



(14) 



and a nonadiabatic coupling H na (t) = -iR _1 (i)^|R(i). If 
the time evolution is slow we can neglect the nonadiabatic 
coupling; then Eq. (fT3"|) is easily solved, 



b a (t) = b a (—oo) exp 



e a {t')dt' (a = +,0,-). 

(15) 

Our major approximation is based on the assumption 
that the population dynamics takes place mainly during 
the time interval when fi -C \JTiT c . i.e. during the 
rising edge of the pump pulse (which is the last to arrive) . 
This assumption derives from our interest in the regime 
of large ionization, which requires strong laser fields. In 
this case we can approximate the eigenvalues and the 
eigenvectors by assuming that the ratio fi / vTTc is a 
small parameter [l6| . After simple algebra we obtain 



£_ 



A + 

n 2 



r 2 - 4A<5 
A + 5- 



n 2 i] + A-5 

4?7 1] + A + T 



n 2 i]-A + 5 

4?7 77 - A - T 



(16a) 
(16b) 
(16c) 



where 
A = 



1-,-%. 



r 



(q + i) v/TT^, 



(17a) 

(17b) 
(17c) 



n 



A-S- 



1 {Yi + r - r c ) 



{q + if Viscid) 



The corresponding eigenvectors are 

n 77 + a - 5 



<p+{t) 

<Po(t) 
<P-{t) 



sin^, cos^, sin^ 
25Q 



nr 



r r) + A + S 

_ n 2 (f 2 + 43 2 ) 

2(f 2 - 4A5) 2 ' f 2 - 4AT f 2 -4A 

n v -A + s n2 

;COs£, — sm4, cos 4 



(18a) 



1 T 

18b) 



Tri-A-5 
where the complex- valued angle £ is defined as 

f 

tan 2^ = ™ 

A-5 



(18c) 



(19) 



Because the pulses are applied in the sequence ionizing- 
control-pump and the population is initially in state rpi, 
the initial adiabatic-state amplitudes are bo(— 00) = 1, 
b±(— 00) = 0. In the adiabatic limit the population 
Po(t) — \bo(t)\ 2 of state ipo(t) evolves as 



Po(t) 







exp 


-iS(t) J' 



n 2 (t' 



-00 r 2 (yj') - 4A(t')5{t>) 
The populations of the original states are 

n 2 {t)[f 2 (t) + AS 2 {t)} 2 



-dt' 



Pi® 

P2(t) 
Pc(t) 



1 - 



2[T 2 (t)-4A(t)S(t)Y 



Po(t), 




Po(t). 



(20) 

P (t), (21a) 
(21b) 
(21c) 



The fluorescence signal is calculated from Eq. ([T]) and 
()21bl) . and then the ionization signal from Eq. ([9]). 

In the following section we will use these formulas to 
examine how to minimize the fluorescence F, i.e. ^2(^)5 
and simultaneously to maximize the ionization /. 



IV. NUMERICAL EXAMPLES 

We compare the analytical results derived in the pre- 
ceding section with numerical simulations for the fluo- 
rescence signal of Eq. ([1}, the ionization signal of Eq. 
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Control Pulse Center (units of T) 

FIG. 5: Contour plots of the ionization signal as a function of 
the control pulse center r c and the control ionization rate r c o. 
The ip 2 loss rate is F = 1/T (top frames), T = 10/T (middle 
frames), or T = 100/T (bottom frames). The Fano parameter 
is q = 1 (left frames), q — 3 (middle frames), or g = 6 (right 
frames). The other parameters are T = T c = T, Tio = 50/T, 
flo = 50/T, ri = — T, 5 = 10/T, A = 0. The number J max 
atop each frame indicates the respective maximal ionization 
signal. Without the control field J max is 0.978, 0.822, and 
0.317 for r = 1, 10, and 100, respectively. 



n=l/T,qBl,I aO.989 r=l/T,q=3,I mas =0.997 r=l/T, q=6, 1 =0.998 




-20 20 -20 20 -20 20 



Control Detuning 8 (units of 1/T) 

FIG. 6: Contour plots of the ionization signal as a function of 
the control laser detuning 8 and the peak control ionization 
rate F c0 . The i> 2 loss rate is T = 1/T (top frames), r = 10/T 
(middle frames), or V = 100/T (bottom frames). The Fano 
parameter is q = 1 (left frames), q — 3 (middle frames), or 
q = 6 (right frames). The other parameters are T = T c = T, 
r ffl = 50/T, n = 50/T, n = -T, t c = -0.5T, A = 0. 
The number l max atop each frame indicates the respective 
maximal ionization signal. Without the control field J ma x is 
0.978, 0.822, and 0.317 for F = 1, 10, and 100, respectively. 



([9]), and the populations (f2~Tj) of states ipi, ip2 and t() c , 
derived from numerical integration of Eq. ([3]) , with the 
Hamiltonian (jHJ. We assume Gaussian pulse shapes, 

fi(t) = n e- [t - T)2 ,T \ (22a) 
I\(t) = T m e-^rt, (22b) 
r c (t) = r c0 e-(*-^) 2 /^ 2 . (22c) 

We use the pump pulse duration T as a time unit and 
1/T as a frequency unit, and choose the center of the 
pump pulse to define the zero reference point of time, 
t = 0. All remaining parameters are variable: the peak 
pump Rabi frequency il , the peak ionization rates T i0 
and r cU , the centers of the ionizing and control pulses Ti 
and t c , their widths Ti and T c , the detunings A and S. 
The Stark shifts are assumed zero because, as we have 
verified, they do not affect significantly the ionization 
signal. For the Fano parameter we have chosen three 
values: q = 1, q = 3, q = 6, which are close to the 
experimental values for LICS in sodium (q — 3.7) [13], 
helium (q = 0.73) [lj|, and hydrogen atoms (q = —5.9) 
Q, for electric-field mixing in rubidium (q = 3.3) 19 1, 
and for configuration mixing in potassium (q = 1) 2C|, 
rubidium (q = 0.1 - 0.3) and cesium (q w 0.43) [2l(. 

Figure [5] shows contour plots of the ionization signal 
as a function of the center of the control pulse r c and the 



peak control ionization rate r cU for different values of the 
Fano parameter q and the irreversible loss rate T from 
state ip2- Larger Fano parameters are clearly favourable 
for ionization but improvement is seen for q = 1 too. In 
principle a lower Fano parameter could be compensated 
by stronger ionizing and control fields. For r cU = (near 
the horizontal axis in each frame), the plain STIRAP 
into continuum, discussed in Sec. UH occurs. The bot- 
tom frames demonstrate that for strong loss (FT ^> 1) 
the LICS-STIRAP improves ionization dramatically, for 
example, from 0.317 without the control field to 0.875 
with it. We have verified that for larger laser intensities, 
a nearly complete ionization can be achieved. 

Figure [S] shows contour plots of the ionization signal 
as a function of the control laser detuning S and the con- 
trol peak ionization rate r cU . Again, the presence of the 
control laser pulse, and the ensuing LICS, are essential in 
achieving high ionization signal, even for strong loss rate 
from state "02 • The asymmetry of the ionization signal 
vs the detuning is typical for the Fano LICS profile. 

In Fig. we display the fluorescence and ionization 
signals vs the timing of the ionizing pulse, for different 
values of the irreversible loss rate T and different widths 
of the control and ionizing pulses. The figure demon- 
strates that efficient ionization requires a counterintuitive 
pulse ordering, with the ionizing laser applied before the 
pump laser (r^ < 0, with optimum about Tj = — T). It 



7 



1.0 



0.6 



\ 1 /> 








- y 


y 

t 

i" 


■ \ 


t 


l 

n 


A ; 


it \ 


l\ i 




'v i i i - 



1.0 - 






0.8 - 


\ / 

V / 




0.6 - 




w 
* 


0.4 - 




* 


0.2 - 




r V i 


- 







1.0 



-6 -4 -2 2 4 6 



-6 -4 -2 2 4 6 



1.0 



0.4 








* 

* A 

* / 

1 


/ F 
t 

- 


i . 


i 

1 i 








-6 -4 -2 2 4 6 



-6 -4 -2 2 4 6 



Ionization Pulse Center t. (units of T) 

FIG. 7: (Color online) Ionization I and florescence F signals 
vs the center n of the ionizing laser pulse. The ip2 loss rate is 
F = 10/T (top frames) or T = 100/T (bottom frames). The 
pulse widths are T* = T c = T (left frames) or Ti = T c = 3T 
(right frames). The other parameters are q — 3, F^o = r c o 
50 Y . Sl = 50/T, t c = -0.5T, A = 8 = 10/T. The solid 
curves show the analytical results derived from Eqs. JT]), ([9]), 
and (I21[) . the dots show the numerical results. 



also shows that an increase of the pulse widths (right 
frames) leads to broadening of the ionization profile but 
does not affect appreciably the maximal ionization signal. 
The figure evidences an excellent agreement between the 
analytical theory and the numerical simulations. This 
agreement indicates that indeed, the ionization dynam- 
ics occurs during the rising edge of the pump pulse, when 
the ionizing and control pulses are already present. The 
manner in which the pulses terminate is not important, 
as evident when comparing the left frames (where the 
pump pulse starts and terminates later) and the right 
frames (where the pump pulse start later but terminates 
earlier because its width is shorter than the others) . This 
is one of the main differences between conventional STI- 
RAP between discrete levels Q and LICS-STIRAP pro- 
posed here. In conventional STIRAP both the initial and 
final times of the pulses are important (the pump must 
start and terminate last). In LICS-STIRAP the initial 
times are important, but the final times are not, because 
there is no population left in the discrete states. In this 
respect, LICS-STIRAP is similar to STIRAP between 
lossy states [Ta | . 

Figure [5] shows the fluorescence and ionization signals 
vs the timing of the control pulse. The figure demon- 
strates that the best timing between the ionizing and 
pump pulses is at r c ~ Tj/2, as used in other figures. 
However, the technique is relatively robust against the 
control timing and moderate deviations from this pre- 
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FIG. 8: (Color online) The ionization and florescence signals 
vs the center r c of the control laser pulse. The loss rate from 
state i> 2 is F = 100/T and the pulse widths are T, = T c = T. 
The other parameters are q = 3, F^o = r c o = 50/T, f2o = 
50/T, n = -T, A = 5 = 10/T. The solid curves show the 
analytical results derived from Eqs. ([9]), and (|21[) . the 

dots indicate the numerical results. 



scription do not affect the ionization efficiency very much. 
The background signal appearing for large deviations of 
t c from this region is produced by c-STIRAP. It serves as 
a reference for the influence of LICS on STIRAP, signifi- 
cant in this case, once again. This figure reveals another 
excellent agreement between analytical theory and nu- 
merical simulations. 

Figure[5]shows the ionization signal as a function of the 
irreversible loss rate T from state tp2- The three curves 
show how the ionization efficiency decreases with T for 
TPI, c-STIRAP (no control laser), and LICS-STIRAP. 
The c-STIRAP efficiency is clearly superior to TPI, and 
LICS-STIRAP adds further considerable improvement 
over c-STIRAP. Note the horizontal logarithmic scale in 
r. Remarkably, LICS-STIRAP maintains high ionization 
efficiency even when the intermediate state ip2 can decay 
hundreds of times during the interaction. Again, the fig- 
ure reveals an excellent agreement between the analytical 
theory and the numerical simulations. 

Figure [TU] shows the ionization signal vs the Fano pa- 
rameter q and the peak control ionization rate. For mod- 
erate values of q (1 < q < 10) nearly complete ionization 
is achieved for sufficiently strong control pulses. When 
q is too small LICS is not sufficiently strong to simulate 
the presence of a bound state and to create a quasi-dark 
state. When q is too large, a very large atomic coher- 
ence is created through the continuum, so that state ip c 
is directly involved in the dynamics, with some atomic 
population transfered to this state. 
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FIG. 9: (Color online) Ionization signal vs the ■02 irreversible 
loss rate V for three different techniques. TPI: coincident 
pulses ("Ti = r = 0), no control pulse (r c o = 0); c-STIRAP: 
counterintuitively delayed pulses (r^ = —T, r — 0), no control 
pulse (r c0 = 0); LICS-STIRAP: counterintuitively delayed 
pulses (n = —T, t = 0), with control pulse (r c o = 50/T, 
t c = — 0.5T). The other parameters are T = T c — T, q — 3, 
Pio = 50/T, Q = 50/T, A = 5 = 10/T. The solid curves 
show the analytical results derived from Eqs. {TJ, (J5J), and 
(|21|) . The dots display numerical results. 
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FIG. 10: Contour plot of the ionization signal as a function 
of the Fano parameter q and the peak control ionization rate 
Pco- The other parameters are T = T c = T, T i0 = 50/T, 
fio = 50/T, T = 100/T, n = — T, r c = -0.5T, 5 = 10/T, 
A = 0. 



V. CONCLUSIONS 

We have demonstrated, using analytical techniques 
and numerical simulations, that STIRAP can be used as 
a tool for efficient ionization of the population of a dis- 
crete state ipi coupled to a continuum via a lossy state i/'o- 
The ionizing laser must precede the pump laser coupling 
states ipi and "02 , and both lasers must be strong enough 
to enforce adiabatic evolution. The ionization probabil- 
ity is further enhanced when a LICS is created into the 
continuum and STIRAP is directed into this LICS be- 
cause then a quasi-dark state composed of state tpi and 
the LICS is created, via which the population flows into 
the continuum. We have shown that almost complete 
ionization can be achieved even when the lifetime of the 
resonantly coupled state ip2 is much shorter than the laser 
interaction duration. We have also shown that in the adi- 
abatic limit, the main population dynamics takes place 
during the rising edge of the pump pulse, when it is lost 
irreversibly either via ionization or fluorescence. Hence in 
LICS-STIRAP only the order in which the pulses arrive 
is important, a result significantly different from conven- 
tional STIRAP between discrete levels, where the order 
of the pulse terminations is also very important. 

Experimental verification of these results will open op- 
portunities for many applications, such as photoioniza- 
tion of ultracold atoms with efficiency close to unity and 
negligible population into intermediate discrete states. 
It should be noticed that LICS-STRAP relies on atomic 
structures in the continuum reached by appropriate laser 
sources. An accurate examination of the atoms listed 
above and their experimental investigation shows that 
sodium is particularly suitable for experimental verifica- 
tion of the proposed technique. 
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